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ABSTRACT 

The  field  of  gaseous  electronics  has  beoome  one  of  burgeoning 
importance  in  recent  years  due  to  the  vast  strides  made  in  the  field  of 
communications  and  other  related  fields.  Theories  abound,  but  few 
impregnable  mathematical  descriptions  of  the  phenomenon  of  the  positive 
column  have  been  advanced.  An  attempt  was  made  to  extend  one  of'  the 
mathematical  solutions  to  a  more  general  case0  The  resulting  equations 
were  programmed  into  a  digital  computer  which  produced  results  in- 
dicating the  efficacy  of  the  program,  but  more  data  will  have  to  be 
digested  in  order  to  establish  the  reliability  of  the  basic  equations. 
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CHAPTER  I 
INTRODUCTION 

This  thesis  evolved  from  the  analysis  of  a  paper  by  Professors  M« 
Satosi  Watanabe  and  llorman  L.  Oleson  of  the  U.  S.  Naval  Postgraduate 
School    in  which  the  phenomenon  of  traveling  density  waves  in  a 
positive  column  was  investigated.   It  seems  appropriate  in  the  beginning 
to  define  exactly  what  is  meant  by  some  of  the  terms  and  expressions  that 
will  be  used  hereafter.  This  is  done  expressly  to  avoid  confusion  over 
semantics  and  not  as  a  gratuitous  collective  insult  to  subsequent  readers. 

In  the  simplest  instance  let  us  assume  that  we  have  a  cylindrical 

glass  tube  filled  with  a  pure  rare  gas  at  a  pressure  of  say,  12  mm  Hg. 

In  terms  of  every  day  experience  this  would  be  considered  a  relatively 

low  pressure;  however,  most  authorities  in  the  field  of  gaseous  conductors 

refer  to  this  as  a  high  pressure.  Sealed  in  either  end  of  our  gas  tube 

let  us  imagine  that  we  have  two  electrodes  which  conform  to  the  shape  of 

the  tube.  Let  each  of  these  electrodes  be  connected  to  an  external  souroe 

of  e.m.f.  through  an  external  resistor.  How,  it  is  well  established  that 

if  either  the  e^m.f.  or  the  resistor  is  varied,  a  glow  discharge  can  be 

induced  to  occur  in  the  tube.  Generally  speaking,  the  glow  discharge  is 

a  steady  state  gaseous  discharge  whioh  appears  at  reduced  pressures  and  is 

characterized  by  relatively  low  current  densities  and  relatively  high 

potential  gradients;  perhaps  it  is  better  to  say  in  this  latter  case, 

relatively  high  values  of  -,  electric  field  to  pressure  ratios,  about 

P 

which  more  will  be  said  later.  This  gaseous  discharge  has  a  typical 
structure  which  is  marked  by  a  rather  steep  potential  gradient  at  the 


cathode  and  operates  basioally  by  electron  liberation  resulting  from 
positive  ion  bombardment  at  the  cathode.   Quite  simply,  this  amounts  to 
nothing  more  than  having  the  gas  act  as  a  oonductor.  However,  the  pheno- 
mena occurring  within  the  gas  while  it  is  acting  as  a  conductor  cannot  be 
so  simply  explained.   In  fact  so  many  phenomena  occur  that  it  is  not 
feasible  to  discuss  all  of  them  within  the  scope  of  this  thesis.  As 
indicated  in  the  opening  statement,  we  shall  confine  our  attention  to  a 
single  phenomenon  in  the  positive  column  of  the  glow  discharge.   In  order 
to  facilitate  our  development  we  shall  make  some  simplifying  assumptions 
about  the  composition  of  our  positive  column.  We  shall  assume  that  it  is 
composed  of  three  fluids: 

a.  A  negative  ion  "gas"  which  we  shall  further  assume  to  be 

composed  entirely  of  eloctrons 

b.  A  positive  ion  "gas"  whioh  we  shall  further  assume  to  be 

composed  either  of  atoms  or  molecules  of  the  gas  with 
which  the  tube  is  filled,  but  each  of  which  carries  a 
single  positive  charge 

c.  The  neutral  particles  of  the  gas  itself— these  particles 

are  by  far  in  the  majority. 
It  is  customary  in  discussing  gaseous  conductors  to  refer  to  ionized 

regions  containing  high  and  approximately  equal  concentrations  of  positive 

(2) 

and  negative  ions  as  the  plasma   .  With  the  proper  environment  there  can 

be  caused  to  appear  in  the  positive  column  alternate  dark  and  light  bands. 
These  bands  may  move  or  they  may  be  stationary.   In  either  case  they  are 
called  striations.   If  they  appear  to  move,  they  are  called  moving 
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striations;  if  they  do  not  move,  they  are  called  standing  striations» 
As  a  parenthetical  matt or  of  interest,  it  seems  well  to  note  that  moving 
striations  aro  not  detect ablo  by  the  unaided  eye* 

The  traveling  density  waves  referred  to  earlier  may  or  may  not  be 
moving  striations.  Certainly  ono  could  be  excused  if  in  the  first  glow 
of  understanding,  he  impetuously  ascribed  the  title  of  moving  striations 
to  these  traveling  density  waves;  however,  sober  and  mature  consideration 
would  make  such  a  conclusion  doubtful.  The  important  point  to  remember 
here  is  that  moving  striations  are  an  observable  phenomena  and  not  an 
ethereal  speculation. 

A  mathematical  approach  by  Watanabe  and  Oleson  has  been  undertaken. 
In  their  approach  they  attempt  to  provide  a  theoretical  basis  for  ex- 
plaining moving  striations.  However,  they  assert  that  certain  of  their 
approximations  may  be  in  error  if  applied  to  actual  striations   •   It  is 
the  intent  hereafter  to  justify  or  refute  these  assumptions.   In  order  to 
do  this  it  is  proposed  to  extend  the  linear  treatment  of  Watanabe  and 

Oleson  to  the  non-linear  case  and  thereby  arrive  at  a  more  general  solution. 
s 

A  search  of  the  literature  was  instituted  to  garner  the  latest,  most  exact 
data  available.  Where  these  data  were  found  to  be  either  doubtful  or  non- 
existent, data  were  taken  from  current  experiments  being  performed  at  the 
U.  S.  Naval  Postgraduate  Schoolo  The  ultimate  aim  was  to  substitute  these 
data  in  the  pertinent  mathematical  equations  and  solve  them  by  means  of  a 
digital  computer.  The  solutions  then  could  be  compared  with  the  results 
obtained  experimentally  and  a  reasonable  conclusion  could  be  reached  as  to 
the  validity  of  the  original  premises,, 
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Sinoe  we  have  assumed  that  we  have  a  fluid  filled  tube,  let  us  estab- 
lish a  coordinate  system  to  describe  the  particles  in  our  fluid.   In  this 
system  the  neutral  particles  are  always  fixed  in  space  and  only  the  positive 
and  negative  particles  move  with  respect  to  our  coordinate  system* 

If  P  is  the  density  of  any  physical  quantity,  Q  the  rate  per  unit 
volume  at  which  the  quantity  is  being  produced,  and  5 V  the  volume  occupied 
by  the  quantity;  then, 

Let  us  define  t»V  as  a  small  volume  of  fluid  and  agree  that  £  V  always 
designates  a  volume  element  that  moves  with  the  fluid  suoh  that  it  contains 
the  same  fluid  particles.   In  general  the  volume,  &V,  will  then  change  with 
time.   If  we  assume  that  SV  is  in  the  form  of  a  rectangular  box,  then  let 
us  say  that  its  dimensions  are  ix,  Jy,  and  izt    such  that; 

JV-  in  iy  U 
How,  the  x-oomponent  of  fluid  velocity,  vx,  may  be  different  for  the  op- 
posing faces  of  our  box  moving  in  the  x-direction0   If  so,  £x  will  change 
with  time  at  a  rate  equal  to  the  difference  between  these  two  velocities; 
thusly,  i       n 

and  by  a  similar  development  we  can  say, 


htr*%     «J    $i*-*i* 


Therefore, 


J  c~\j 
Inserting  this   in  the  equation  above  for    — —     we  have 


syd?+?v^SV-Q£V 


dt 


dividing  out  *JV, 


jt^V.7,Q 


but by  definition, 


'^^~r7-ve 


therefore, 

Which  by  a  relation  in  vector  analysis  reduces  to: 

This  of  course  is  nothing  more  than  the  equation  of  continuity.  A  more 
rigorous  derivation  can  usually  be  found  in  most  texts  on  mechanics. 


See  for  instance  K.  R.  Symon,  Mechanics  (Add i son-Wesley  Publishing 
Company,  Inc.,  Cambridge,  Massachusetts,  1953),  p.  265  ff» 


For  the  sake  of  simplicity  we  shall  assume  that  all  of  the  negative 
ions  are  electrons  and  that  the  positive  ions  are  formed  by  the  oollision 

of  eleotrons  with  neutral  particles  and  by  no  other  process— a  reasonable 

(3) 
assumption  if  we  know  that  our  gas  is  not  an  electro-negative  gas   ' .   See 

fig.  1  for  a  simple  graphical  representation  of  this  ionization  phenomenon. 

6—  O  ©  ^ 

F'g    I     (Icrsjj&a     Gam,    vov  fv&-*i) 
If  we  place  a  charged  particle  in  an  electric  field,  the  partiole  will 

move  under  the  influence  of  the  Coulomb  law  of  attraction.   It  seems 
reasonable  that  if  there  are  neutral  particles  in  the  field  which  are  un- 
affected by  Coulomb  forces,  the  charged  particles  moving  in  the  field  will 
collide  with  the  neutral  ones.  This  obviously  ooours.  However,  since  the 
charged  particle  is  still  under  the  influence  of  the  field,  the  particle 
will  again  tend  to  move  after  collision  toward  the  electrode  opposite  in 
sign  to  the  partiole.  The  overall  effect  is  to  give  the  charged  particle 
a  motion  generally  in  one  direction  (see  fig.  2). 

Since  the  Coulomb  forces  operate  to  move  the  particle,  it  would 
appear  that  the  same  law  acting  between  oppositely  charged  particles  would 
cause  them  to  recombine.  This  does  happen  under  proper  conditions;  i.e., 
high  ^pressures  (of  the  order  of  one  atmosphere).  However,  at  low  pressures 
the  collision  frequency  and  collision  losses  are-  low  and  the  relative 
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velocity  of  the  charged  particles  is  higher.   Because  the  velocities  are 
higher,  collision  without  recombination  takes  place  since  the  attracting 
field  is  not  powerful  enough  to  deflect  the  ions  and  make  them  move 
closely  together  over  a  large  enough  distance^  '. 
Let  us  now  define : 

n"  s  negative  ion  density 
n'  -   positive  ion  density 
v  =  negative  ion  drift  velocity 
tt  =  positive  ion  drift  velocity 

z  =  number  of  ion  pairs  created  per  unit  time  per  one 
electron  (constant  if  the  temperature  variations 
are  very  small   ') 
then  from  the  equation  of  continuity: 

^  t 
6qn(1)  ur*7-(«'7-)-arfS0 

Further  consideration  of  a  charged  particle  in  an  electric  field 
will  inevitably  lead  to  the  conclusion  that  the  velocity  of  the  ion  is  a 
function  of  the  field*  Experimental  evidence  shows  that  the  velocity  is 
proportional  to  the  field  to  a  first  approximation,  viz: 

(see  figs.  3  and  4) 
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Figure   3,      Drift  velocity  of  electrons   in  argon  as  a  function  of  E/p. 
R.  A.   Nielsen,   Phys.  Rev.   50,    950   (1936) 
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Figure  4«     The   drift  velocity  of  Argon  atomic   ions   in  Argon  as  a  function 
of  E/p   .      J.  A.  Ilornbeck,    Phys.   Rev.   84,    615    (1951) 


The  use  of  vector  notation  is  justified  when  consideration  is  given  to  the 
fact  that  the  motion  of  the  ion  is  directed.  We  may  therefore  say: 

where  K  is  a  proportionality  constant  called  the  mobility — a  quantity  which 
will  be  more  fully  discussed  latere  Considering  again  figs0  3  and  4,  it 
oan  be  seen  that  the  magnitude  of  the  velocity  of  the  electron  is  much 
greater  than  that  of  the  positive  ion;  indeed,  in  certain  instanoes  the 
velocity  of  the  electron  can  be  shown  to  be  a  thousand  times  as  great  as 
that  of  the  positive  ion.  This  fascinating  facet  of  the  properties  will 
be  given  more  light  later  in  the  disoussion  of  ambipolar  diffusion* 

Envisioning  a  cloud  of  ions  (or  electrons)  in  the  tube,  one  is  led  to 
the  inevitable  conclusion  that  the  charged  particles  may  exist  in  varying 
densities •   If  this  be  the  case,  then  obviously  the  ions  will  flow  from 
regions  of  higher  density  to  the  regions  of  lower  density  by  diffusion. 
This  will  of  oourse  contribute  a  component  to  the  velocity  at  which  the  ions 

move,  and  we  oan  say: 

— »• 

where  7*1  is  the  concentration  gradient  with  the  negative  sign  indicating 
that  the  motion  occurs  in  the  direction  of  decreasing  concentration. 
Therefore, 

7  -  -  ^  vn 

where  D  is  a  proportionality  constant  called  the  diffusion  coefficient  for 
the  gas,  A  more  adequate  treatment  of  the  derivation  of  this  last 
equation  can  be  found  in  some  detail  in  most  texts  dealing  with  the  subject 
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of  diffusion*   Combining  these  contributions  to  velocity  we  have, 

eqn  (2)   ?=KF-|^ 

The  electric  field  under  which  the  particles  act  is  composed  of  two 
parts,  the  impressed  eleotric  field  and  the  electric  field  generated  by 
the  space  charges;  such  that, 

E'-E.+  E' 

where  E0  is  the  impressed  eleotric  field  and  will'  be  considered  constant, 
and  E*  is  the  electric  field  generated  by  the  space  charges  and  depends  on 
the  number  of  electrons  and  ions  present;  i.e.,  the  density.  Therefore 
we  can  say  (from  Poisson's  Equation): 

V'E  -?.(£.  +  £')  =  v-E'*fire(n+-n) 

where  e  is  the  electronic  charge. 

It  is  proposed  to  solve  equations  (l)  by  assuming  that  the  actual 
solution  is  a  superposition  of  a  steady  uniform  solution  n*  plus  a  non- 
steady,  nonuniform  disturbance  jf»      In  order  to  linearize  the  equation 
V/atanable  and  Oleson  assumed  that  the  nonsteady  nonuniform  disturbance  was 
small  compared  to  the  uniform,  steady  disturbance^,  which  is  not  in  agree- 
ment with  experimental  results.^  '  ^    '      In  the  following  development  no 
such  assumption  will  be  made— -to  the  oontrary,  we  shall  assume  that  at 
some  particular  timp  the  orders  of  magnitude  of  the  steady,  uniform  and 


See  for  instance  L.  3.  Loeb,  Basic  Processes  or  Gaseous  Electronics 
(University  of  California  Press,  Berkeley  and  Los  Angeles,  1955),  p. 188  ff. 
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nonsteady,  nonuniform  parts  may  be  approximately  equal.  This  is  not  to 
imply  that  they  are  always  equal,  but  rather  that  they  could  be  equal  at 
some  particular  time.  The  word  "uniform"  as  used  here  means  that  the 
quantities  do  not  depend  on  the  x-coordinato  taken  along  the  axis  of  the 
tube^'.   In  keeping  with  our  proposed  solutions  let  us  now  write: 

r\r-  r\l  +  jf 

Qqn  (3)    .. 

n  -  n6tp 

We  define  £  as  the  field  generated  by  the  steady,  uniform  charge 
densities  n£. 

Substituting  the  value  given  for  V    in  equation  (2)  into  equation  (l) 

^  +  V-  irf  (-Kt  -  {%.  V "})]  -  zri-o 

at. 


£+  7-  [^(k'B  '  {-£  vn>})] - ZrT*D 


clearing  the  parentheses: 
at  l— 

^!t  v>  [n+Kr£-(Dfvnr}]-itn--=o 

Substituting  equation  (3)  in  the  preceding  equations: 

Directing  our  attention  to  'the  terms  containing  n*  and  recalling  that  by 
definition  these  are  constant;  therefore, 
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which   leaves, 
—  +  v  [-n'tK£  -  7j'  K'£  -  D~  V  v]  -  2  n;  -  2z>'?  o 

d  t 

dt  L- 

but  we  have   stipulated  that  E  a  E0  •/    E1;    it  then  follows: 

|^  +V'[;rilKE0-  n; lis'- u'K~E0-v~l< 'z'-D' VvJ-zn;-? zT=a 

Again  we   refer  to  the  fact  that  E0  is   a  constant,   consequently. 
Omitting  the   zero  terras,    our  equations  become: 

It  r  ?'  C  ^  K  £'-#  *"£0-  ^~  tffe'  -  Vvd]-  2ft„-2  jfzo 

3  t 

£t 

We  have  defined  £  as  the  field  generated  by  the  uniform,  steady  charge 
density  nQ.   Let  us  now  define  ~e   as  the  field  due  to  the  nonuniform, 
nonsteady  disturbance  jf  ;  such  that, 

t--£te 

Inserting  thi3   relation  into  the   equations   above: 

j^  +v-[:n\Ki- nl  \(^- tj' KE6-  v K'i- v K  e-  if  v p'2' ^ril- 2v  -o 
2t        L 
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Removing  brackets, 

*v  -  „  -  *   - 

^*i;K7-£-i.K7'?-^£4 P-zT-  KvV-E-  K£-  Vtj'-Ktj  v-  e  - Ke-  Vv 

From  Poiss on's  Equation: 

where  e  is  the  charge  on  the  electron.  Recalling  that  our  region  of 
interest  is  near  the  axis  of  the  tube  removed  from  the  ends  and  assuming 
cylindrical  symmetry,  we  find: 

2* 


Yfith  these  conditions   in  mind  let  us   rewrite  the  equations: 

7i 


- Kif  [h ire  (z/-  jJfr-KZ  |£  - 0'  p{-  £r?; -  2 v  =o 


tf&birelTf-vfl+K'rtW-jfW.^^o 


Inasmuch  as  our  interest  in  tho  field  is  in  the  central  region  along  the 
axis  of  the  tube,  we  note  that  the  field  £  (or  more  particularly  6^  )  will 
be  zero  because  of  the  steady,  uniformity  of  the  charge  density  n*  causing 
it.  Further,  in  this  region  n*  is  approximately  equal  to  n". 

ii -  J)'  pi- K£0  |£"  +  If, V-  Lpr-K~ez &-W-v)v- zK=o 


|£-d*  ^V^v^-i^S  ^'4^y-^ 
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where, 

)2~- trre K~ n;  ^-HireX 

Ij-tircK'Wt)  t=  Hrrel<f 

For  comparison  purposes  the  terms, 

which  appear  in  equations  (4)  do  not  appear  in  the  original  paper  by 
.Vatanabe  and  Oleson.  These,  then,  are  the  terms  dropped  by  vYatanabe  and 
Oleson  as  a  result  of  the  mathematical  expediency  employed  by  them  to 
linearize  their  equations,  and  as  a  result  of  disregarding  certain  higher 

order  terms. ^  ' 

Equations  (4)  are  the  ones  for  whioh  constants  will  be  computed  for 

a  particular  set  of  conditions;  e.g.,  Argon  gas,  12  mm  Hg  pressure,  eto« 
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CHAPTER   II 
EVALUATION  OF  CONSTANTS 
Reiterating  our   original  contention,  we  propose  to  apply  equations    (4) 
to  an  experimental   situation.      In  order  to  do  this  certain  constant  terms 
must  be   evaluated.      Other  constants  are  well  known  but  will  be   repeated 
here   in  order  that   our  results  may  be   identically  reproduced.      In  this 
latter  category  of  constants  we  find: 
e   =  4.8  X  10"10  e.s.u. 
k  =   1.38  X  10"16  ergs/  °K 
tf  =   3.1416 
In  the  former  class   of  constants  there  will  be  a  need  to  know  the  values   of 
K'#  K",    D7,    D  ,   E0,    z,   n*,    and  nj.     Rather  than  list  them  arbitrarily,    it 
seems  more   appropriate  to  derive  them  individually  in  order  to  indicate  their 
reliability. 

The  particular  experimental   situation  to  which  we   shall  apply  equations 
(4)   is   as  follows: 

A  tube   0.5  cm  in  radius  containing  pure  Argon  at   a 

pressure   of  12  mm  Kg  with  currents   of  from  20ma  to 

( G ) 
lOGma  in  our  gaseous  conductor^    ' „ 

With  these  conditions  the  value   of  the   applied  electric   field  oan  be   found. 

From  fig.   5     it  can  be   seen  that   for  a  pressure-radius  product   of  6  we  have 

in  our  case   an  -     ratio  of  approximately  0.5.     At   a  pressure   of  12  mm  Hg 

this  will  yield  an  EQ  of  6  volts/cm.     The  value   of  EQ  which  we   obtained  for 

our  case  yields  a  value   of  2  X  10       statvolts  when  converted  to  e.s.u. 


For  a  description  of  the  method  of  obtaining  these  curves  refer  to  the 
source  cited  in  connection  with  figure  5. 
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From  probe  measuromonts  in  the  plasma  within  the  tube  it  has  been 
experimentally  determined  that  ^°/, 


P 


10 


O.I 


- .  J . 


ao 


nQ  =  n'Q  i   7  X  10xu  ions/cm" 


2  5^'d 

200/v 

o  z  4  & 

Figure  5 
(for.i?ed  Gases ,  von  kngel) 


£ 


/£> 


It  appears  rather  obvious  that  a  relation  exists  between  mobility, 
diffusion,  and  the  number  of  ion  pairs  created  per  unit  time  per  one 
electron.   Since  the  latter  two  terms  have  been  defined,  it  remains  only 
to  define  the  first  term — mobility.  V/e  have  rather  arbitrarily  stated 
that  mobility  is  the  proportionality  constant  which  relates  the  drift 
velocity  to  the  field.  However,  it  seems  more  important  to  have  a 
"feeling"  for  the  phenomenon  of  mobility  than  to  have  a  pedantic  knowledge 
of  the  word  itself.  Consider  a^ain  a  charged  particle  moving  under  the 
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influence  of  a  field  and  within  a  region  oooupied  by  the  neutral  partioles 

of  aA(fig.  2).  We  oan  say,  then,  that  mobility  is  the  drift  velocity  of 

* 

an  ion  moving  in  a  unit  field. 

In  order  to  utilize  the  experimental  data  available^6',  it  is  neces- 
sary to  determine  values  of  mobility  for  the  positive  ions  and  the 
electrons.   In  dealing  with  positive  ion  mobility  it  was  found  that  there 
is  a  difference  in  mobilities  depending  on  whether  the  ions  present  are 
"atomic"  ions  or  "molecular"  ions*   It  has  been  shown  that  when  pressures 
of  the  order  of  6  mm  Hg  are  encountered  in  Argon,  A£  and  A'   ions  are  found 

* 
In  searohing  for  a  best  value  for  mobility  the  phrase  "reduced 
mobility"  was  often  encountered.  Reduced  mobility  is  best  described 
as(7); 

■where  p  is  the  gas  pressure  and  T  is  the  gas  temperature  in  °K,  Many 
results  are  given  in  terms 
fig,  4),  This  reduced  pre: 


results  are  given  in  terms  of  reduoed  pressure,  p  ,  (see  for  example 

1  pressure  is  defined  as  l®): 


273p 

Po  =  -y* 

Very  often  in  studying  graphically  presented  results  one  observes 
that  experimenters  and  theoreticians  will  plot  ion  mobility  versus 
field  to  pressure  ratio   (i.e.,  §)  —where  the  pressure  may  be  reduoed 
pressure  or  that  of  the  gas;       P  e.g.,   figs.  3  and  4.      It  oan  be 
shown(2)  that  p  is   inversely  proportional  to  the  mean  free  path, 
thusly:  ( 

P~X 

If  this  be  so,  then: 

■where  EA     is   an  expression  for  the  energy  acquired  in  a  m.f.p. 
This  gives  the  researoher  a  more  vivid  and  accurate  picture  of  what 
is  happening. 
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to  be  present   in  the  ratio  of  0.5^    '.     The  mobilities  of  these  two  ions 

wore  found  to  be   2.6  and  1.6  ?Wgfig       respective ly(9).     Assuming  that  the 

volt /cm  r  J  o 

two  types   of  ions   are  distributed  in  the  ratio  indicated  throughout  the 
positive  oolumn,    a  weighted  average  mobility  was  found  to  be: 

KT=   4.04  X  104  om2/stv-seo 
If  the  two  types  of  ions  are  assumed  to  be  present   in  equal  concent rat ions, 
then: 

tT  =  4.38  X  104  cm2/stv-seo 
As  oan  be   seen,  the  orders  of  magnitude  are  the  same  or  at  least  near 
enough  that  either  torm  could  be  used  with  some  degree  of  accuracy. 

Prom  fig.  3  the  drift  velocity  of  an  electron  was  determined  to  be 
0.36  X  10     cm/sec.     To  a  first  approximation, 

r-  -  v  -  0.36  X  106   .  ,  0  Y  ,n7   2/  . 

K  -  -  -  ■  ■■  '   -  1.8  X  10  cm  /stv-seo 

2     2  X  10-2 
The  diffusion  coefficients,  D'  and  D~,  are  most  easily  found  by 
employing  the  mobility.  The  relationship  between  mobility  and  diffusion 
is  derived  in  numerous  texts  on  the  subject^2'.  The  relation  is: 

at       -i  ~ 

D       kT 
where  o   is  the  electronic  charge,  k  is  Boltsaaann's  constant,   and  T  is  the 
ion  temperature   in  °K.     Tho  temperature  of  the  ions  is  approximately  equal 
to  the  temperature  of  tho  gas^      '•     From  this  equation,  then,  the 
following  values  wore  determined: 

DTs  3.48  om2/seo 

D"  S  6.46  X  104  cm2/sec 
In  order  to  compute  the  value  of  z  it  is  necessary  to  employ  ambi- 
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polar  diffusion  in  tho  positive  column  results  from  the  previously  shown 
fact  that  the  diffusion  coefficient  for  electrons  is  greater  than  for  ions. 
Assume  that  as  an  original  condition  there  exists  near  the  axis  of  the 
positive  column  an  almost  equal  number  of   electrons  and  positive  ions. 
After  a  short  interval  the  eleotrons  will  tend  to  diffuse  toward  the  walls 
of  the  tube  at  a  rate  greater  than  the  positive  ions.  This  will  cause  a 
negative  charge  to  be  set  up  on  the  walls  of  the  tube.  The  preponderance 
of  positive  ions  remaining  near  the  axis  of  the  tube  (due  to  their  slower 
diffusion)  causes  an  electric  field  to  be  established  between  the  axis  of 
the  tube  and  the  wall.  This  field  in  turn  will  cause  a  slowing  down  of  the 
olectrons  moving  toward  the  wall  and  a  speeding  up  of  the  ions.   It  happens 
that  wall  recombination  is  preponderant  at  the  pressure  we  have  specified. 
This  is  borne  out  by  the  absence  of  recombination  spectra  indicating  that 
volume  recombination  is  not  important  here.   In  an  equilibrium  state  the 
losses  due  to  recombination  are  equal  to  the  rates  of  production.  The 
electric  field  across  the  positive  column  is  such  that  electrons  arriving 
from  the  faraday  dark  space  are  accelerated  to  produce  enough  electrons 

and  ions  per  unit  length  of  the  tube  to  compensate  for  losses,  thus 

(2) 

establishing  equilibrium^  '.  Specifically,  in  our  case  we  have  restricted 

our  consideration  to  the  region  near  the  axis  of  the  tube  so  that  radial 
variations  need  not  be  considered.  The  point  to  be  remembered  here  is 
that  under  the  action  of  the  radial  electric  field,  an  average  velocity 
of  diffusion  is  soon  reached  such  that, 

4  dn"   dn^ 

V   —  V   ■  V         .  1   -  — — 

r  "  r    rav»   dr   '  dr 
which  leads  to  the  fact  that  tho  total  ambipolar  current  is  zero.  As  we 
have  already  seen, 
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Also,  wo  have  previously  implied  (page  13  )  that  near  the  axis  of  the 

tube  e"  ■  0  and  6X  -  0,     Using  the  expressions  for  velocity  and  eliminating 

the  field  between  them,  it  can  be  shown  that^  ': 


D^"*  -/  D'TST 


a  =  T^  K" 

where  D&  is  defined  as  the   ambipolar  diffusion  coefficient.     Further, 
since  K"   is  much  greater  than  K  '  and  T      is  muoh  greater  than  T      (where  T( 
and  Tp  are  the  electron  and  ion  temperatures  respectively),    it  can  be 
shown  that : 


=   DT    o  ■  D     KT-KTJ 


Da;DT_e.D     IT   -  KT e 

Tp  ir  e 

Employing  known  data    '6)    (Te  s   12,550  °K  and  Tp  s  300  °K)   in  the   approxi- 
mate formula  for  Da  yields: 

Da  =  145  cm  /sec 
From  the   ambipolar   solution  to  our  problem  it  can  be   shown  that: 

*          (2.405)2  n 
z   = u 

(R)2        a 

where  R  is  the  tube  radius o   Inserting  the  value  for  D&  computed  above 

( 6 ) 
and  the  radius  of  our  tube— 0»5  cm  '  for  the  conditions  we  have  chosen, 

we  find: 

z  -   3,35  X  10  ion  pairs/electron-sec 

Inserting  the  constants   just  computed  into  equations    (4)  yields: 

dt  djt  £-/.  r  <};* 
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Little  accuracy  is  sacrificed  by  disregarding  those  values  of  constants 
■which  are  dominated  by  the  exceedingly  largo  values  of  other  constants 
which  are  additive „ 

In  order  to  solvo  this  rather  formidable  differential  equation, 
little  error  is  involved  if  the  partial  derivatives  are  approximated  by 
total  derivatives  and  minutely  small  increments  are  taken  to  be  used  in 
a  digital  oomputer.  This  may  be  more  clearly  seen  later  where  the  above 
equations  are  set  up  for  programming  into  the  computer.   Once  the  equations 
are  written  in  a  "useful  form"  for  the  computer,  we  find  we  must  have  some 
starting  point;  i.e.,  some  initial  values  for  v^,   x,  ex.,  etc.  at  some 
time,  t.  Since  the  approach  we  have  taken  appears  to  be  a  rather  unique 
one,  there  were  no  readily  available  precedents  by  which  we  could  be 
guided.  Harking  back  to  the  original  intent  of  this  thesis,  the  answer 
suggested  itself.   Since  V/atanabe  and  Oleson^  '   assumed  a  solution  to 
their  linear  equation,  it  seemed  reasonable  to  take  their  evaluation  and 
apply  it  to  our  problem.   If  this  satisfied  our  equations,  then  credence 
is  lent  to  the  original  hypothesis  thereby  extending  the  theoretical  con- 

« 

tention  to  actual  conditions.     Watanabe   and  Oleson      '   assumed  linear 
solutions: 

V    n   a/eikx-itJt 

X>*=   a""eikx"i<jJt 
where   aT  and  a"   are  complex  amplitudes,   k  is  the  wave  number,    and  cj   is  the 
angular  frequency.     Thus,    another  problem  arose — the  feeding  of  a  complex 
number   into  a  digital  computer.      It   is  not   impossible  to  execute   suoh  a 
program,   but    in  view  of  the  physical  construction  of  the  computer,    it  would 
be  necessary  to   sacrifice   some   of  the   accuracy  of  the   original  program  in 
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ordor  to  do  so.  However,  since  our  solution  consists  of  a  real  part  and 
an  imaginary  part;  that  is  to  say,  a  real  solution  and  an  imaginary 
solution,  either  is  a  solution  to  our  problem.   So,  we  will  program  only 
the  real  solution.  This  reduces  our  solution  to: 

if  -  a~cos  (kx-ot ) 

2/ r   a7c  o s  (loc-6Jt ) 
Letting  the  time  arbitrarily  be   zero   initially  these  reduce  to: 
iJ0  m  a~cos  kx0 
7Jf «  a/cos  kx0 
Again  we  have   introduced  terms  which  must  be  evaluated.      Directing  our 

attention  to  k,   the  wave  number,    and  referring  to  the  work  of  Watanabe 

(l) 
and  01esonv    ',  we   observe  that  we  must   assume  their   linear  counterparts 

of  our  equations    (4).     This  results   in  no   inconsistencies  when  the 

restrictions   placed  on  our  equations   are  considered.     Watanabe   and  Oleson's 

equations   are:  ^ 

whore  K,  Y*   /j  and  ^  are  "the  same  as  the  ones  previously  defined. 
Substituting 

V  z   a-eikx-^ 

*  -  4   ikx-O; 

V  -   a'e 

into  Watanabe   and  Oleson's  equations   and  performing  tho    indicated 
operations  will   give: 

(-iw-ih-Eok  ■/■    D'k2  ■/  I  )a~  -  Hz  or/  =   0 
-X,a"  4  (-1^7^    iK/E0k  /    D/k2  ■/   13)aJ  =   0 
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In  order  for  Wo  homogeneous  equations  to  have   a  solution  the   determinant 
of  coefficients  must  be   zero* 


=   0 


-i<U-iK"E0k  ■/  D"k2  -/i  -  y2 

Upon  solving  this   determinant  we  have, 

[-iW-   iK"EQk  ^D~k2  -/     £]    [_-i6J  yiK^E0l:  7^    d42  Y3]  -   Y*  ^   =   0 
Recalling  that   if  a  complex  equation  can  be  written: 

R  ■/■    il  =  0 
where  R  is  the  real  part   and  il  the   imaginary  part,   then  R  -  0  and  il 
for  the  equation  to  be   satisfiedo 

More  particularly  for  our  case, 
for  Rao, 

-O2  ^K74l0k-UK"E0k-/K"K^E2^:27^  (D"k2-/    tfj  )  -  *-  )fv        =0 
for   1=0   (since   i  /  0), 

-w[D/kV(j7  -K~E0k[Vk2  /K3]    -o)[D"k2  ^tfj  -/K^0k|Vk2  /JfJ 
Solving  this   last   equatidn  foroO: 


=   0 


a  0 


6J  = 


kEQ  {Jc2(k/d"  -  K"D"^     -    (K"  Jj  -  K7^  Y,    )] 
[V    y    D")k2/    (    tf(    /    ^)] 
Eliminating   OJ   between  the  two  equations  yields, 

eqn   (6)  (D"k2  /    V,    XD^k2    ■/   Jfj   )  -  Uy 

/  k2E§(K^    V    K-)2(A2  ^    U  )(D"k2   •/  *,) 
[(d42      /     ^  )   ^     (D-k2   •/     K,   )]  2 

Introducing  new  variables, 

eqn    (7)        ^  =  D"k2  -/  ^      and    ^    =   d42  7^    ^ 
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=   0 


2 

Solving  for  k  in  the  first  of  equations  (7), 


k2  = 


.  s'-v, 


Substituting  this  value  for  k  in  the  second  equation  (7), 

Talcing  cognizance  of  the  fact  that  D"  ^>>  D'  ,  we  can  say: 
Referring  to  tho  original  expression  for  n   , 

and  considering  the  preceding  approximation,  we  observe: 

d42  k   0 

Since  J)T   is  not  zero  (see  page  18  ),  this  approximation  can  be  true  if  and 
only  if, 

k2  t   0 
Inserting  the  approximate  values  of  ri     and  £    from  equation  (7),  et  seq, 
into  equation  (6)  and  observing  that  k  *  0, 

If  we  map  equation  (6)  in  the  ?-r>  plane,  we  have  approximated  very  closely 

a  point  on  the  extremely  complicated  curve  represented  by  equation  (6); 

2  (1) 

i.So,  the  point  where  k  s  0.   It  can  be  shown    from  an  elementary 

geometrical  relationship  that  knowing  a  point  and  the  slope  of  a  curvo  at 
the  point,  the  equation  of  the  tangent  can  be  written.   In  our  case, 


-(8)     r-i-[j$zss]($-8W 
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where,    gs  £&  and  h  =  ^LzLipM 

Substituting  the  values   of   n     and    €    as   shewn  in  equations    (7)  for  the 
appropriate  terras   in  equation   (8),    and  noting  that: 

z    <^    4  7TeK'n 

2 

we   solve   for  k    :  , 

,2  -     (K^     ■/    K")z 

k       -     ■*■-"/"    '*  '    ; j*    q 

47Tne 
Recalling  the  proposed  linear   solution,  we   direct   our  attention  to  the 
complex  amplitudes   and  write: 

&f    =  Ae"1^     (with  a  and  0  real) 


which  can  be  written: 

a7~    -  A(oos  0  -   isin  0) 


but   as  has  been  shown, 

(-i  U)  -  iK"EQk  4  D"k2  7^   ^   )a~  -    \z  a/     -  0 
From  these   last  two  equations  we  can  say, 

J-       »          </        -a    •      «f         D"k2    /    r,     -   i(u)-/    K"i:0l:    ) 
ii     -  Ac os  j2  -   lAam  jo  3     L-~ ■ * J 

a"  K2 

Equating  reals   and  imaginaries, 

rf  -   D"*2   7^  IT.        and     Asin  0  -  "  V    ^Bol: 
Ac  os  p  =  ' r  y 


y, 


Solving  for  A, 


a  s  D-yz  y  1       p-k2  /  y,    jTT^y 
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From  which, 

,       CJ  -f    K"Eok 
tan  p  =   — 

D-k2    4  *, 
Inserting, 


(V  /  lT)k2   /(    Jf,     V  ^)] 


and 


Simplifying  yields, 

,       (JST    J    K~)kEo 

eqn     9         tan  0  =  1 — t-^ — 5—- = 

[(d-/  4  D-)k27/(  jf,  4  y3)] 

It  oan  be   shown (k-L'  that    ban  $  has   one  maximum  between  0  and    y2# 

Talcing  the   derivative   of  tan  p  with  respect  to  k  and  equating  to   zero  we 

find  at  this  point: 


eqn   (10)     k6  = 


2  _  W. 


[2Eo(K"/    /    K")j     (D^   7^    D-)2 
Squaring  equation   (9)   and  inserting  equation   (10)   gives: 

(tan^)2 ^    ^  K')2S° 

w  "  [icn/  /  d-)(  y,  y  y,  )J 

Calculating  this  value   of  ft  from  the  known  values   of  the  constants 
in  the  preceding  equation,  we  used  this   as  the  maximum  value  of  a  phase 
factor  with  which  to  enter  the  computer  along  with  our  previously  desoribed 
boundary  conditions.     The  value   of  $  obtained  as   described  was  found  to 
be   1°54'. 
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CHAPTER  III 

ASSUMPTIONS  AND  CR  IT  I  QUE  OF  ASSUMPTIONS- 
SOME  DIFFICULTIES  ENCOUNTERED 

The  experimental  data  yield  a  Langnuir  plot  which  is  a  straight  line 

( 6 ) 
over  two  oyclos  leading  to  an  assumed  Llaxwellian  distribution  for  electrons v  I 

The  assumption  that  the  ions  have  a  llaxwellian  distribution  stems  from  the 

faot  that  at  several  millimeters  pressure,  the  temperature  of  the  ions  is 

approximately  equal  to  that  of  the  gas  molecules^   '•  The  assumption  that 

we  have  radial  ambipolar  diffusion  under  our  stated  conditions  within  the 

tube  is  based  on  the  rosults  of  Karge,  Hooks,  and  Oleson^  '.   In  this  work 

they  state  that  the  instantaneous  field  intensity  is  due  to  a  very  small 

surplus  of  positive  ions  or  electrons.  Their  reasoning  for  this  is  very 

amply  stated  in  their  work.  Since  we  can  say  that  n"  A  n'  ,  it  follows 

that  we  can  say  the  ambipolar  current  at  the  walls  is  approximately  zoro. 

Further, 

dn"    dn '      -    -4 
-rr—  z     -rr-    and  v  =  v' 
dt     dt       r    r 

whero  we  define  v~  and  v£  as  the  negative  (i.e.,  electron)  ion  and 

t 

positive   ion  radial  drift  velocities  respectively. 

In  our  preceding  work  and  the   subsequent  computer  work  it  was  assumed 

that   z  was  a  constant.     If  z  is  also  evaluated  from  the  relation^    ': 

£  ■  /kT  ^  2 
eqn   (11)      z  *   9  X  107  ape"  lt\^l     V± 

(6) 

using  the  temperatures  observed  by  Karge,  Hooks,  and  Oleson   ,  it  is 

found  that  z  varies  widely  with  temperature  as  shown  below: 

Te 

16,440 

14,660 

8,200 

27 


1, 

,89 

X 

105 

5, 

,3 

X 

104 

2, 

►  22 

Farris    '  notes  a  variation  in  2  based  on  temperatures  (the  data 
being  taken  from  Pupp»s  work).  When  this  information  is  applied  by  the 
uso  of  our  equation  (11 ),  the  following  results  are  obtained: 

Te 

28,000  2.7  X  107 

11,520  1.63  X  103 

It  should  be  noted  that  the  validity  of  the  equation  for  2  depends  on 

having  Ilaxwellian  distribution  of  elootrons  -which  we  have  assumed. 

Farris     has  recently  published  a  short  paper  on  a  theory  of  moving 

striations  based  on  the  diffusion  theory.  His  basic  assumption  is  quoted : 

If  it  is  now  assumed  that  the  small  difference  between  the 
ion  and  electron  concentrations  has  the  same  effect  axially 
as  radially,  and  that  moan  values  of  mobilities  and  diffusion 
coefficients  may  be  used,  the  theory  of  a  stationary  positive 
column  may  be  readily  generalized  to  include  changes  of 
concentration  axially  and  with  time.  Eliminating  the  ourrent 
and  field  from  the  usual  conservation  and  current  density 
equation  gives:  S  -  dn/dt  -  Da  ^n,  vdiere  S  is  the  rate  of 
production  of  ions  per  unit  volume  per  second  and  Da  the 
Schottky  ambipolar  diffusion  coefficient  given  by: 
Da  =  kT€K//e. 

The  major  assumption  different  from  our  assumption  is  that  the  axial 

drift  velocities  of  the  electrons  and  positive  ions  are  assumed  equal 

within  the  positive  column.  Using  our  equations  (l)  and  oldjainating  the 

field  wo  obtain:  nv  =  Da^7n.  Substituting  this  in  our  equation  (l)  we  get 

Farris'  equation  (l).  However,  we  question  the  validity  of  the  assumption 

that  the  axial  drift  velocities  of  all  ions  in  the  positive  column  are 

equal.  The  interesting  point  is  the  variation  of  2  with  electron 

temperature o 

Frara  these  considerations^- °'   wo  observe  that  there  is  a  fallacy  in 

our  assumption  that  z  is  a  constant.  True,  we  imposed  the  condition  that 
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this  was  valid  only  for  small  temperature  variations.   But  at  the  time 
we  imagined  this  variation  oould  easily  be  on  the  order  of  100  °I1.  How, 
it  appears  that  this  is  not  so.  Certainly  this  is  a  field  for  further 
inquiry.  How  this  deviation  from  our  original  assumption  affects  our 
problem  will  be  shown  in  the  next  chapter.  We  have  taken  zn^  to  be  a 
constant  of  the  value  of  2,345  X  1014, 

It  seems  noteworthy  to  point  out  that  there  is  a  phase  difference 
between  the  maximum  degree  of  ionization  and  the  maximum  degree  of 
excitation  which  has  been  observed  in  some  rare  gases,  however  this  was 
not  observed  by  Karge,  Hooks,  and  Oleson*  ',  As  a  point  in  fact,  they 
observed  that  a  maximum  in  light  intensity  was  accompanied  by  a  maximum 
TG,  and  a  minimum  in  light  intensity  was  accompanied  by  a  minimum  Te. 
This  would  suggest  that  there  is  no  phase  difference  in  the  case  of  Argon 
ions  in  Argon  gas  at  12  ram  Hg  pressure.  If  we  can  assume  that  oquation  (11 ) 
holds  and  further  assume  that  the  maximum  light  intensity  is  also  evidence 
for  maximum  excitation,  then  we  can  say  indeed  there  is  no  phase  differ- 
ence between  the  phenomena. 

The  temperature  dependence  of  mobilities  and  diffusion  coefficients 
was  ignored  in  this  thesis.   Since  the  evaluated  constants  in  equation  (4) 
were  taken  to  be  average  values^  ',  it  is  suspected  that  temperature 
dependence  is  of  little  significance. 

The  relationship  between  mobility  and  diffusion  used  in  this  thesis, 

K    e  ,  (page  18)  may  be  a  tenuous  one.  Von  Engel^  '  states  that  the 
D  "  kT" 


* 

Private  communication  with  II.  L,  Oleson,  who  with  A,  V/,  Cooper,  has 

found  that  maximum  light  intensity  and  maximum  electron  tenporature  are 
not  coincident  in  neon  at  several  mm  prossuro. 
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the  equation  is  not  applicable  to  electrons  for  at  least  two  reasons; 

electrons  drift  with  velocities  which  are  not  proportional  to  the  applied 

field  and  their  mean  free  path  depends  on  their  speed*  In  defense  of  our 

use  of  this  equation  we  reliod  on  experimental  data  (see  fig.  3)  which 

shows  that  for  our  conditions  the  velocity  is  approximately  proportional 

to  the  appliod  field.  The  E     obtained  by  Karge,  Hooks,  and  Oleson  was 

"p  max 
6  volts  per  cm  per  mm  IIg,   This  figure  is  well  within  the  range  of  fields 

wherein  experimental  data  reveals  that  velocity  is  proportional  to  field 

(fig,  3),  No  attempt  has  been  made  to  justify  or  refute  the  m,f,p, 

(a) 

difficulty.  Von  Engelv  '  further  states  that  the  equation  holds  in  the 

oase  of  ions  which  move  and  diffuse  in  moderately  strong  fields:   in 
strong  fields  the  velocity  is  proportional  to  the  square  root  of  the  field, 
He  further  states  that  it  can  be  assumed  that  the  equation  holds  for  weak 
fields.  That  our  conditions  fall  within  the  purview  of  Von  Engel's  state- 
ments can  be  seen  from  figure  4,  wherein  the  velocity  of  the  ions  is  shown 
to  be  proportional  to  the  field, 

Graham  and  Ruhlig^12'  obtained  a  value  of  2  X  10  cm  /sec  for  the 
diffusion  coefficient  of  the  elootron  (D")  by  using  Droyvesteyn  distri- 
bution in  conjunction  with  the  Boltzmann  transport  equation.  This  valuo 
was  calculated  at  the  same  conditions  specified  for  our  case.  As  can  be 
seen  from  page  18  we  oaloulated  a  value  of  6,46  X  104  cm  /sec.  No 
attempt  was  mado  to  resolve  this  difference.  However,  as  will  be  shown 
later,  this  difference  may  bo  rather  easily  investigated  by   use  of  the 
program  we  propose  by  merely  varying  the  coefficients  of  certain  terms. 
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CHAPTER   IV 

PROGRAMMING 
In  order  to  facilitate   an  understanding  of  the  computer  program  we 
will   set   down  a  brief  mathematical  description  of  what   is  taking  place 
within  the  computer.      For  convenience   let    V    =   z   and    V  =  w—  equations    (4) 
may  then  be  written: 

o  X 

•/  ^w  /  2.345  X  1014 
ii.      lf=   D"ji|y:K-(E0/ex)i^.   [I   -     4TTeK-(Z-w)]    w 


■/  /  z  V  2.345  X  10 


14 


and  from  Poisson's  2qn, 

iii.     < — ^  s  4lfe(z  -  w) 

dX  ' 

Furthermore,  we  arbitrarily  establish  as  an  initial  condition: 

ex  -  0,  at  x  e  0 

Let  us  now  divide  the  part  of  our  positive  column  wherein  our  equations 

will  apply  into  n  equal  parts. 

Let : 

x  denote  tho  origin  of  our  column  (  =  0) 

X-,  be  the  first  point  at  a  distance  Ax  from  xQ 

then, 


Let: 


x^  -   i£x,  where   i  s   0,    1,    ,   n. 

and, 

x  a  x^  -  xi-l  a  Wn»    L  being  the   length  of  the   column 

Zi    3     Z(xi,t) 
Wi    a    W(xi,t) 

ex-s  ©x(xi)  ■  4  7re(zi  -  w^)        (from  iii.) 

31 


Remembering  that  it  was  proposed  to  replace  the  partial  derivatives  with 
total  derivatives  (page  21  )»  w<3  now  observe  the  technique  of  doing  this 
by  the  finite  approximation  formulae: 


u.  /-   -  u. 


au,     s  .  1/1    i 

\2  u.  /.  —  2u.  J.  u.  , 
^  x2  (Ax)2 

Utilizing  this   expediency  on  our  first  three  equations   in  this  chapter, 
we  have : 

i,a)         d?  =  (aT^Vi"  2zi  ^Zi-l}  "  Ix^o^XiXVi-  Zi) 


n*a 


-["!(,  7^4ireK^(Zi  -  w^Z-l  /    ^Wi/2.345  X  10 

dw_. 


14 


dw.        n-  „- 

)  _i  -  !> -(w,   /.-   2w,    «/w.    -)  t^— (En  ■/    ey    ft     ,    -  w.  ) 

dt  (ax)2     ^-T1  x         1"1         Ax     °  xi^i^l        i 


-[y,   -   4fTeK~(z.    -  w.)]wi   •/     ^Zi  -/ 2.345  X 


10 


14 


iii.a)  ev     =  ev         -/4ire(Ax)(s*   -  w- ) 

xi         xi-;L  1  1 

It   is   proposed  to  examine   our   solution  over  ten  increments;   therefore, 
we  restriot   our  problem  as   follows: 

15  1,    2,    3,    ....    ,    10 

wQ=   zQ  -  0 

wn=   zn  =   0 

Rearranging  equations  i.a),  ii.a),  and  iii.a)  above,  we  have: 


">     ^^-^0^^)]^, 
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>l4 


{&o  4^)  "  g^  "  |V   4^(zi  "  *i>]}  2i 


fit* 


+  [^rj  zi-l  4h^4  2.345  X  io14 


ii.b) 


dt" 


2L-  /E(e_  /ex   )lwi  y. 


■/r~~J^i-l  4    \H  4  2.345  X  1014 


iii.b) 


Jx4  : 


_(ax)2J 
eXi_1  4  4^e(Ax)(2i  -  Wi) 


Let  us  now  introduce  some  new  constants  by  combining  same  of  those  whioh 
frequently  oocur;   specifically  let. 


ai  = 


s  s   .387     (arbitrarily  setting  x  a  3cm) 
A?  =  EL  a  1.35  X  104 


(ax)' 
K7< 


to      ■»  ' 

z  -  4  x 


kV  : 


A3  s  4TteK~    =   2.44  X  10 


-4 


Bl  = 


(AX) 


2     s   .720  X  104 


K" 


B2  -  ~-  =   .6  X  101 

6     *"      A    X     "" 


A  X 


B3  =  4ffeK"  -  .109 
Insorting  these  new  constants   into  equations  i.b),   ii.b),   and  iii.b) 
yields: 

1,0  }       d?*["Al  "  A2(E°  ^i^i/l 

/^A2(E0  -/eXjL)  -   2AX  -  A3j^;  •/    (Ei  -  w^JJ     2; 
7^1zi-l  4   h^i  4  2.345  X  1014 
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ii,c)    at*  a  [Bl  "  B2(E°  ^  ^i^i/l 

/  Blwi-1  -/B3nozi  V    2«345  x  lo14 

iii.c)     eY     =  ex     ,   -/  6.023  X  10~9(ax)(z.   -  w.) 
xi         xi-l  x         x 

It  should  be  noted  that  the  value  of  n~  is  arbitrarily  takon  to  be  7  x  10  0 

Yfhile  this  may  or  may  not  be  exact,  it  is  the  correct  order  of  magnitude. 

Further,  in  the  ambipolar  region  n"  a  n£  •   In  ordor  to  avoid  overcomplicating 

the  manipulation  of  decimals  a  new  variable  is  introduced: 

VT±    -    1015w! 

Changing  variables  and  inserting  the  values  of  the  constants  where  applicable: 


+  .327  ^'■_)  t.noi  *l0*\j'.  +  ,23i5 


+ .  720  X  IV  Wjm  +  .  743  X/o'V-  f-,  2JV5 


e*r  ev/  ♦"■'w^c*;-^') 


34 


Let  us  define  s-jj   -   .02  -f    ex     and  So.   =  z-   -  w.,  then  the  preceding 
equations  oan  be  written  as: 

^1  ~-[.W-,\35  *I0S  StiJgl+{l35  M<fs,r  ,771  -.?mrf*[.7  W+  s^jl; 
+.38 1  i!     f  ,1769  X/O*  ty  +  ,Z3¥S 


'0      ' 


+  , 726  t  tor*;    +.7tf*'6  2;t,23i5 


H 


e^  -eximl  +  .i807xio*s2j 


These  are  the  final  forms  of  the  equations  which  the  computer  will  be 
asked  to  solve.  A  general  idea  of  the  method  of  solution  can  be  obtained 
by  considering  these  equations  in  conjunction  with  the  flow  diagram  which 
appears  on  the  following  page.  The  details  of  execution  of  the  program 
are  found  following  the  flow  diagram. 
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yes 


•*i 


0 


,  1807  t  lOtSz\ 

fit,-.,  t,/gQ7//Q*32l-  -e>,' 


oz  *■  e^i-  s, 


./35  X/05  5u' 


.389-  ./35  X/o5Si{ 


a(--D89-./35//055/i)^',>, 


■  7/tO^Szi 


-.244  KlO'2(.7/ilO",t  Szi) 


-JVi  -[.244XiO*(  7)<t0'*t  S2i)] 


.'35  X/05S,j  -{.77H  +  [244  *lO,*(.7)UO"'t  Sjii)]) 


Sv  -  ,1708  X  I09vJ- 


a.ta* 


a3t(a,'  32) 


a*  *£>,■»■(  a, -r  a*  )J 


2*i-z,23«s  +  [aitCa3*(al  +  a2j)} 


.  6  x  /o7  S// 


.720%  10*-  .LX'O7  $ti 


6,  -  C  7*0  X  10"  ~.LXi09S,i  )*',-. 


■  7x/o'"-S2i 


J09x/o'5(.ix/&"'-S2;) 


I.Hi  XiOH  +  [109  xiO's(.7xtO'"-  S*/)J 


,LXlO*S,it{l.H'miO',t[lOiXIO'S(.7*iOm',-Ssi)] 


^h/tXEtX/Q7S,;»O.V»X/D^Qo^/fl,5(.7X/0-,,-Sa/)]}3 


-«'— I 


b3-.720*IO*  X  ^j-i 


fcv  =  .743  X/0/0   ?J 


4.*L 


*>j+(6^b,) 


bH+[b3  +  (b.+bt)} 


%j-=..23*s+{b,+[b3  +  a>li-b*)]} 


H-*J 


rest    »>?/•/    >  J 


tfO 


Nutfietitil    Integration     Sob-  routine    (Rvnge-  Kutta-  Gift  method) 
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The  following  is  the  detailed  program  of  the  computer  solution  to  our 
problem. 
Loo  Opn 


Opn 
Code 


M 


0010  md 

0011  ad 

0012  ad 

0013  ad 

0014  su 

0015  tm 

0016  md 

0017  tm 

0020  sf 

0021  mr 

0022  ad 

0023  ad 

0024  tm 

0025  ad 

0026  ad 

0027  tm 

0030  ad 

0031  ad 

0032  tm 

0033  md 

0034  tm 

0035  ad 

0036  mr 


26  2100 

35  2100 

35  0200 

35  0210 

36  I"  z± 
34  0202 
26  2100 

34  3000 
31  0240 

25  0241 

35  2100 
35  2100 

34  3000 

35  2100 
35  2100 

34  3000 

35  2100 
35  2100 
34  2001 

26  2100 

34  3000 

35  2000 
25  2001 


1,1 
2 

2100 

0201 

0200 

0207 
i 

Wjj 

2100 
2100 
2100 
0202 
2001 
0205 
0206 
2100 
0242 
0243 
2100 
2000 
2001 
2100 
2100 
2100 
0244 
0245 


M 
3 

0205 

0200 

0207 

0014 

O202J 

0020 

2000 

0022 

2000^) 

200li 

2002") 

2003J 

0400 

2002^\ 

2003  J 

0400 

0203  -\ 

0204J 

0035 

2000 

0037 

2000  "^ 

2001 J 


Remarks 


0-*e 


xi 


1  X  2-12  /  1  X  2-24  /  i  X  r12-/  i  X  2'24 


zi  "  wi  B  s2i  in  °202 

no -> 001 6,  yes-»002CT 

n  (2000 
u   (2001 


-.0™  v  m8  v      •   (2000 
.1807  X  10°  X  sti  in  (2001 


e     in  (2002 
xi-l    (2003 


gives  ex.  in  [2gg0 

02  in  S2002 
•u«j  in  j2003 


(2000 
xi  =  *li  ^   (2001 


.02  «/e 


sn  in  0203  and  0204 

no  -*0033,  yes -"0035 
0  —-  2000,  2001 


-.135  X  10  sli 
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0037  ad 

0040  ad 

0041  tm 

0042  tm 

0043  md 

0044  tm 

0045  ex 

0046  ad 

0047  ex 

0050  tm 

0051  md 

0052  tm 

0053  tm 

0054  ad 

0055  ad 

0056  ad 

0057  ad 

0060  ad 

0061  ad 

0062  tm 

0063  tm 

0064  md 

0065  tm 

0066  ad 

0067  mr 

0070  ad 

0071  ad 


35  2100  0246  2002 

35  2100  0247  2003 

34  3000  2100  0400 

34  2001  2100  0045 

26  2100  2100  2000 

34  3000  2100  0054 
32  0014  0211  0212 

35  0212  0213  2004 
32  2004  0211  0050 
34  [zi/l]  2100  0053 
26  2100  2100  2000 
34  3000  2100  0054 

34  3000  2100  0500 

35  2100  2000  0214 
35  2100  2001  0215 
35  2100  2100  2000^ 
35  2100  0202  2001^ 
35  2100  0250  2002" 
35  2100  0251  2003J 
34  3000  2100  0400 
34  2000  2100  0066 
26  2100  2100  2000 

34  3000  2100  0072 

35  0252  2000  2000' 
25  0253  2001  2001. 
35  2100  0254  2002> 
35  2100  0255  2003> 


.389  in  2002  and  2003 

.389  -  .135  X  105  sn  in  [200I 

no  -*0043,  yes-0045 


no-*0051,  yes -0053 


a,  i-  (2000 

al  ^  (2001 

ai  in  °214 

1  (0215 


sai  in  2000  and  2001 


.7  X  10"4  in  2002  and  2003 

.7  X  10"4  /  sai  in  2000  and  2001 
no-- 0064,  yes -0066 


.244  X  1012  (.7  X  10"4  *f  s£i) 


.774  in  2002  and  2003 
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0072  tm 

0073  ex 

0074  tm 

0075  md 

0076  tm 

0077  tm 
0100 

0101  ad 

0102  ad 

0103  ad 

0104  ad 

0105  su 

0106  ex 

0107  tm 

0110  md~ 

0111  tm 

0112  tm 
0113 

0114  ad 

0115  ad 

0116  ad 

0117  ex 

0120  tm 

0121  md 

0122  tm 

0123  tm 


34  3000  2100  0400 

32  0212  0211  0074 

34  [  z±  1  2100  0077 

26  2100  2100  2000 

34  3000  2100  0101 

34  3000  2100  0510 

35  2100  2000  0216^ 
35  2100  2001  0217> 
35  2100  0257  2001 

35  2100  0212  2004 

36  2004  0213  2004 
32  2004  0211  0107 

34  [si-l]  210°  °112 

26  2100  2100  2000 

34  3000  2100  0114 

34  3000  2100  0516 

35  2100  2100  0220 
35  2100  2001  0221 
35  0212  0226  2004 
32  2004  0211  0120 
34  "  T  wj_  1  2100  0123 
26  2100  2100  2000 
34  3000  2100  0125 
34  3000  2100  0524 


-.774  -.244  X  1012  (.7  X  10"4  /  s£i)  in 
2000  and  2001 


no -*-0075,  yes-*  0077 


,   (0216 
a2  **   (0217 


.387  in  2001 


no-*0110,  yes-*' 0112 


a3  in  0220  and  0221 
a3  in  0220  and  0221 


no -»  0121,  yes  ^0123 
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0124 

0125  ad 

0126  ad 

0127  ad 

0130  ad 

0131  tm 

0132  ad 

0133  ad 

0134  tm 

0135  ad 

0136  ad 

0137  tm 

0140  ad 

0141  ad 

0142  tm 

0143  ad 

0144  tm 

0145  ta 

0146  pr 

0147  ht 

0150  sm 

0151  sm 

0152  ex 

0153  ad 

0154  ad 

0155  tm 


35 
35 
35 
35 
34 
35 
35 
34 
35 
35 
34 
35 
35 
34 
35 
34 
33 
21 
22 
30 
30 
32 
35 
35 
34 


2100 
2100 
2100 
2100 
3000 
2100 
2100 
3000 
2100 
2100 
3000 
2100 
2100 
3000 
2000 
3000 
0157 
2000 

f 
2001 
0014 
2001 
0154 
2002 
3000 


2000 
2001 
0220 
0221 
2100 
0216 
0217 
2100 
0214 
0215 
2100 
0262 
02G3 
2100 
0227 
2100 
2000 
0147 

f 
2000 
0227 
0156 
0277 
2100 
2100 


0222' 

0223. 

2002' 

2003- 

0400 

2002 

2003 

0400 

2002" 

2003, 

0400 

2002' 

2003. 

0400 

2000 

0145 

0150 

0002 

f 
2002 
2004 
0154 
0154 

f 
0160 


a4  in  0222,  0223 


a3  in  2002,  2003 

a3  /  a4  in  2000,  2001 


a2  in  2002,  2003 


a  2  */  a3  /  a4  ^   2000,  2001 


*l   in  2002,  2003 


a1/.a2/a3  /a4  in  2000,  2001 


.2345  in  2002,  2003 


£fi  in  2000,  2001 

^i  (^1  scaled)  ^  2000,  2001 


dz-i 

— t  xn  2002 

df 

[  zi]  inm3  of  2004 


fzH  in  m3  of  0122 
dzj. 


3T~ 


transferred  to   integration  sub-routine 
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0156 

f 

f 

7777 

0157 

f 

f 

0001 

0160 

inn 

34 

0204 

2100 

0163 

0161 

md 

26 

2100 

2100 

2000 

0162 

tm 

34 

3000 

2100 

0170 

0163 

ad 

35 

0203 

0264 

2000 

0164 

mr 

25 

0204 

0265 

2001 

0165 

su 

36 

2100 

2001 

2001. 

0166 

ad 

35 

2100 

0266 

2002 

0167 

ad 

35 

2100 

0267 

2003 

0170 

tm 

34 

3000 

2100 

0400 

0171 

tm 

34 

2001 

2100 

0174 

0172 

md 

26 

2100 

2100 

2000 

0173 

tan 

34 

3000 

2100 

0311 

0174 

ad 

35 

0120 

0213 

2004 

0175 

ex 

32 

2004 

0211 

0304 

0176 

"tan 

34 

3000 

2100 

0304 

0177 

0200 

0201 

0001 

0001 

f 

0202 

0203 

0204 

0205 

0206 

0207 

no-»-0161,  yes-*0163 


.6  X  10 '  s1JL  in  2000,  2001 


•720  X  104  in  2002,  2003 


.720  X  104  -.6  X  10'  sii  in  2000,  2001 
no  -> 0172,  yes  -►0174 


} 


i  X  2"12  /  i  X  2 


s2j_  temp  storage 


!li 


-24 


2  X  (i  X  2"12  /  i  X  2"*24) 


41 


0210  su 

0211 

0212  mr 

0213 

0214 

0215 

0216 

0217 

0220 

0221 

0222 

0223 

0224 

0226 

0226 

0227 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0240 

0241 


36   1305  1306  0202 

7777  f  f 

25     f  2001  2001 

0002  f  f 


0001 
02     f 


dummy  for  0014 
extracted  for  0034 
dummy  for  0036 

an  temporary  storage 
&2  temporary  storage 
a3  temporary  storage 
a^.  temporary  storage 


dzi  temporary  storage 
3F" 


f    f 
f   0030 


f   0033 


1047     3476     0000 


b"L  storage 


bg  storage 


b3  storage 


b/   storage 


.1307  X  108 


42 


0242 
0243 
0244 
0245 
0246 
0247 
0250 
0251 
0252 
0253 
0254 
0255 
0256 
0257 
0260 
0261 
0262 
0263 
0264 
0265 
0266 
0267 
0270 
0271 
0272 
0273 


02 


02 


02 


f 
0121 

f 
3227 

f 
3071 

f 
0000 

f 

3441 

f 
6142 

f 
3061 

f 
1011 

f 
1700 

f 
2670 
0000 
1604 

f 
3062 

f 
3407 


f 
7300 

f 
3777 

f 
2601 

f 
2226 

f 
3567 

f 
2335 

f 
1156 

f 
1730 

X 

4061 

f 
6600 
0000 
0000 

f 
1212 

f 
7777 


f 

0000 
0017 
7720 

f 
0142 

f 
3173 
0047 
6333 

f 
1362 

f 
4571 
0033 
0000 

f 
1156 
0030 
0000 
0017 
0000 
0060 
2176 
0017 
4512 


.02 


-.135  X  10' 


.389 


.7  X  10~4 


-2.44  X  10 


-.774 


12 


.387 


G 


.1708  X  10 


.2345 


.6  X  10' 


.720  X  104 


.109  X  1015 


1.44  X  104 


43 


0274 

0275 

0276 

0277 

0300 

0301 

0302 

0303 

0304 

tan 

34 

0305 

md 

26 

0306 

tan. 

34 

0307 

tm 

34 

0310 

0311 

ad 

35 

0312 

ad 

35 

0313 

su 

36 

0314 

ad 

35 

0315 

ad 

35 

0316 

ad 

35 

0317 

tia 

34 

0320 

tm 

34 

0321 

md 

26 

0322 

tin. 

34 

0323 

ad 

35 

0324 

mr 

25 

0325 

ad 

35 

0326 

ad 

35 

f  f  0041 

7066  2107  7440. 

0013  0013  f 

f  f  0027 


fwi/ll  210°  °307 

26   2100  2100  2000 

34   3000  2100  0311 

34   3000  2100  0540 


2100  2000  0230 

2100  2001  0231 

2100  0202  2001 

2100  2100  2000 


.763  X  1010 

m  -/l  m/l  f,  m  =  12  facta}) 
(aa^l)  X  2"36 


no  -*0305,  yes-*0307 


bx  in  0230,  0231 


sai  in  2001 


.7  X  10 


-4 


s2i  in  2000,  2001 


no  -*321,  yes  -*323 


2100  0250  2002") 

>        .7  X  10"4  in  2002,  2003 

2100  0251  2003  J 

3000  2100  0400 

2001  2100  0323 

2100  2100  2000 

3000  2100  0330 

2000  0270  2000' 

2001  0271  2001. 
2100  0272  2002^) 
2100  0273  2003  J 


as 


-4 


.109  X  10lb  X  (.7  X  10"*  -  s2i) 


1044  X  104  in  2002,  2003 


44 


0327 

tm 

34 

3000 

2100 

0400 

0330 

tm 

34 

0204 

2100 

0333 

0331 

md 

26 

2100 

2100 

2002 

0332 

tan 

34 

3000 

2100 

0335 

0333 

ad 

35 

0203 

0264 

2002 

0334 

mr 

25 

0204 

0265 

2003 

0335 

tm 

34 

3000 

2100 

0400 

0336 

ex 

32 

0120 

0211 

0337 

0337 

"km 

34 

[wi] 

2100 

0342 

0340 

md 

26 

2100 

2100 

2000  A 

0341 

tm 

34 

3000 

2100 

0344  / 

0342 

tm 

34 

3000 

2100 

0550  y 

0343 

f 

0344 

ad 

35 

2100 

2000 

0232 

0345 

ad 

35 

2100 

2001 

0233  j 

0346 

ad 

35 

2100 

0266 

2000^ 

0347 

ad 

35 

2100 

0267 

2001 J 

0350 

su 

36 

0120 

0213 

2004 

0551 

ex 

32 

2004 

0211 

0352 

0352 

tm 

34 

[»I-i] 

2100 

0355 

0353 

md 

26 

2100 

2100 

2000 

0354 

tm 

34 

3000 

2100 

0357 

0355 

tm 

34 

3000 

2100 

0560 

0356 

0357 

ad 

35 

2100 

2000 

0234^ 

0360 

ad 

35 

2100 

2001 

0235  J 

1,44  X  loV-  .109  X  1015  X  (.7  X  lO""4-^) 
in  2000,  20001 
no-»-0331,  yes -0333 


no —0340,  yos-»0342 


b2  in  0232,  0233 


♦720  X  104  in  2000,  2001 


no  -0353,  yes-*0355 


b3  in  0234,  0235 


45 


0361 

ad 

35 

2100 

0274 

2000^ 

0362 

ad 

35 

2100 

0275 

2001  J 

0363 

ex 

32 

0074 

0211 

0364 

0364 

tm 

34 

[4] 

1  2100 

0367 

0365 

md 

26 

2100 

2100 

2000 

0366 

tin. 

34 

3000 

2100 

0371 

0367 

tm 

34 

3000 

2100 

0570 

0370 

0371 

ad 

35 

2100 

2000 

0236") 
0237  J 

0372 

ad 

35 

2100 

2001 

0373 

ad 

35 

2100 

0234 

2002^ 

0374 

ad 

35 

2100 

0235 

2003  J 

0375 

"fcm 

34 

3000 

2100 

0400 

0376 

tm 

34 

3000 

2100 

0440 

0377 

• 

0400 

si 

27 

3000 

0422 

2004 

0401 

ad 

35 

0423 

2004 

0410 

0402 

tin 

34 

3000 

2100 

0426 

0403 

su 

36 

2002 

2000 

2004 

0404 

cm 

30 

2003 

2004 

2003 

0405 

ad 

35 

2003 

2001 

2001 

0406 

to 

37 

2001 

3000 

0415 

0407 

sf 

31 

2000 

2001 

2000 

0410 

0411 

su 

36 

2000 

2002 

2004 

0412 

sm 

30 

2001 

2004 

2001 

.763  X  1010  in  2000,  2001 


no-*0365,  yes  -0367 


b4  in  0236,  0237 


D3  in  2002,  2003 


D3  ■/  b4  in  2000,  2001 


exit 
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0413 

ad 

35 

2002 

2100 

2000 

0414 

tm 

34 

3000 

2100 

0405 

0415 

si 

27 

2001 

0424 

2001 

0416 

sm 

30 

2001 

0424 

2001 

0417 

si 

27 

2001 

0425 

2001 

0420 

ad 

35 

0425 

2000 

2000 

0421 

tm 

34 

3000 

2100 

0407 

0422 

02 

f 

f 

0030 

0423 

■fern 

34 

3000 

2100 

f 

0424 

02  - 

f 

f 

0001 

0425 

f 

f 

0001 

0426 

"tan 

34 

2001 

2100 

0431 

0427 

tm 

34 

2003 

2100 

0433 

0430 

tm 

34 

3000 

2100 

0410 

0431 

sf 

31 

2000 

2001 

2000 

0432 

"tm 

34 

3000 

2100 

0427 

0433 

sf 

31 

2002 

2003 

2002 

0434 

ta 

33 

2002 

2000 

0411 

0435 

tm 

34 

3000 

2100 

0403 

0436 

0437 

0440 

ad 

35 

2100 

0232 

2002 

0441 

ad 

35 

2100 

0233 

2003 

0442 

tm 

34 

3000 

2100 

0400 

0443 

ad 

35 

2100 

0230 

2002 

0444 

ad 

35 

2100 

0231 

2003, 

b2  in  2002,    2003 


b2  •/  b3  /b4  in  2000,   2001 


bi  in  2002,   2003 
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0445  tm  34 

0446  ad  35 

0447  ad  35 
0450  tm  34 


0451 
0452 
0453 
0454 
0455 
0456 
0457 
0460 
0461 
0462 
0463 
0464 
0465 
0500 
0501 
0502 
0503 
0504 
0505 
0506 
0507 
0510 


ad  35 

ta  33 

pr  21 

ht  22 

sm  30 

sm  30 

ex  32 

ad  35 

ad  35 

ad  35 

tm  34 

tm  34 

tm  34 

sm  30 

ox  32 

sf  31 

ad  35 

mr  25 

tm  34 
02 

sm  30 


3000  2100  0400 

2100  0262  2002' 

2100  0263  2003 . 

3000  2100  0400 


2000 
0157 
2000 

f 
2001 
0337 
2004 
0461 
2002 


0227 
2000 
0454 

f 
2000 
0227 
0156 
0277 
2100 


2000 
0455 
0002 

f 
2002 
2004 
0461 
0461 

f 


0201  0200  0200 

0276  0200  0012 

3000  2100  1100 

3000  2100  0010 

0050  0506  2005 

2005  0507  0502 

2100  fzi/.]!  2002 

2000  2002  2000" 

2001  2003  2001. 
3000  2100  0054 

f  f  ,  0014 

f  7777  f 

0212  0506  2005 


bi  /■   b2  -/  b3  •/■  b4  in  2000,  2001 


•2345  in  2002,  2003 


■g—     in  2000,  2001 

dwj_ 
dT~ 


dw-i 
.,„   transferred  to  integration  sub-routine 

(i^l)X  2-l2^(i/DX  2-24^1X2-1V 
JLX2-24 

test  m  ■/  1  y  i 

U.T.  to  integration  sub-routine 


[zi/ll  in  ^2  of  2005 


&1   in  2000,  2001 


48 


0511  ox  .  32  2005  0507  0512 

0512  sf  31  2100  [  ij.1  20°2 

0513  ad  35  2000  2002  2000 ^ 

0514  mr  25  2001  2003  2001 J 

0515  tm  34  3000  2100  0101 

0516  sm  30  0107  0506  2005 

0517  ex  32  2005  0507  0520 

0520  sf  31  2100  fli-lj  2002 

0521  ad  35  2000  2002  2000^ 

0522  mr  25  2001  2003  2001  J 

0523  tm  34  3000  2100  0114 

0524  sm  30  0120  0506  2005 

0525  ex  32  2005  0507  0526 

0526  sf  31  2100  Fv[J  2002 

0527  ad  35  2000  2002  2000^ 


a2  in  2000,  2001 


a3  in  2000,  2001 


a4  in  2000,  2001 


0530  mr   25   2001  2003  2001  J 

0531  tm   34   3000  2100  0125 

(cells  0532  through  0537  are  filled  with  zeros) 

0540  sm   30   0304  0506  2005 

0541  ex   32   2005  0507  0542 

0542  sf   31   2100  W[^i]   2002 


0543  ad   35   2000  2002  2000" 

0544  mr   25   2001  2003  2001 

0545  tm   34   3000  2100  0311 
0546 

0547 


bx  in  2000,  2001 


49 


0550 

sra 

30 

0337 

0506 

2005 

0551 

ex 

32 

2005 

0507 

0552 

0552 

sf 

31 

2100 

[wi] 

2002 

0553 

ad 

35 

2000 

2002 

2000' 

0554 

mr 

25 

2001 

2003 

2001 

0555 

tm 

34 

3000 

2100 

0344 

0556 

0557 

0560 

sm 

30 

0352 

0506 

2005 

0561 

ex 

32 

2005 

0507 

0562 

0562 

sf 

31 

2100 

[wi-l] 

2002 

0563 

ad 

35 

2000 

2002 

2000 

0564 

mr 

25 

2001 

2003 

2001 

0565 

"fern 

34 

3000 

2100 

0357 

0566 

0567 

0570 

sm 

30 

0364 

0506 

2005 

0571 

ex 

32 

2005 

0507 

0572 

0572 

sf 

31 

2100 

[h] 

2002 

0573 

ad 

35 

2000 

2002 

2000 

0574 

mr 

25 

2001 

2003 

2001 

0575 

"bm 

34 

3000 

2100 

0371 

0576 

0577 

b2  in  2000,  2001 


b3  in  2000,  2001 


50 


1000  si  27  3000  1016  2001 

1001  ad  35  1014  2001  1013 

1002  ad  35  2000  2100  1022 

1003  bl  05  3000  3000  1014 

1004  ex  32  2002  1000  2000 

1005  mr  25  2002  2007  2002 

1006  ex  32  2002  2001  2000 

1007  ex  32  2100  2001  2002 

1010  si  27  2001  2005  2001 

1011  tm  34  2001  2100  1005 
1212  pr  21  2000  1020  0001 


1013 

1014 

tm 

34 

3000 

2100 

f 

1015 

02 

f 

f 

0004 

1016 

02 

f 

f 

0030 

1017 

5000 

f 

f 

1020 

01 

f 

f 

f 

1021 

7400 

f 

f 

(cells  1022  through  1047  are  filled  with  zeros) 

1050  ad   35   2100   2100   1057     clear  i  tally 

1051  ad   35   1060   1057   1052 


1052  [ad   35   1305 /i  2100   2000J 


1053  tm  34  3000  2100  1000 

1054  ad  35  1062  1057  1057 

1055  tm  34  1061  1057  1051 

1056  tm  34  3000  2100  1160 


U.T.  to  convert  and  print  sub-routine 
i  -f\  -*i 


51 


1057 

1060 

ad 

35 

1305 

2100 

2000 

1061 

0030 

,  f 

f 

1062 

0001 

f 

f 

i  X  2"12 


m 


■/l    ff 


(The  following  is  the  sub-routine  for  the  integration  of  n  1st  order 
equations  by  the  Runge-Kutta-Gill  method) 


1100  si 

1101  ad 

1102  tm 

1103  tm 

1104  ad 

1105  ad 

1106  ad 

1107  ad 

1110  ad 

1111  [sm 

1112  ad 

1113  j^su 

1114  ad 

1115  [mr 

1116  su 

1117  to 

1120  mr 

1121  sm 

1122  ad 

1123  ad 


27 

35 
34 
34 
35 
35 
35 
35 
35 
30 
35 
36 
35 
25 
36 
37 
25 
30 
35 
35 


3000 
1246 
1163 
3000 
2100 
2100 
2100 
1250 
1234 


1232 

*  1 
1267 

2004 

2004 

2006 

2006 

1251 

1251 

1272 


1242 
2000 
1303 
2100 
2100 
2100 
1177 
1251 
1252 

BJ 
1251 
2001 
1250 
AO 
2005 
2100 
1301 
1242 
1270 
1271 


2000  ~| 

1160  -I 

1200 

1107 

1247 

1250 

1251 

1252 

1111 

200lJ 

1113 

20047 

1115 

2005] 

2006 

1226 

2007 

1270 

1271 

1124 


Prepare  exit  to  main  program 

Test  if  first  time  through 

U.  T. 

0  -»k   (clear  print  tally) 

0  -#►  j    (clear  iteration  tally) 

1  X  2  (reset  var.  no.  tally) 
i  X  2-12  ^    jx  2"24  tally 

[9 J  ]    in  Mx  and  [Vjin  M,  of  1111 

ik\]  in  Mi  of  1113 
[A']    in  M2  of  1115 

test      /^'/>/      (truncation  test) 

i  X  2"36  from  i  X  2"12 
i  X  2"12  7    IX  2"36 
fx\  J    in  t^  and  M3  of  1124 
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1124  fad 

1125  sm 

1126  ad 


1130  fsu 

1131  ad 

1132  fmr 

1133  sm 

1134  sm 

1135  ad 

1136  mr 

1137  ad 

1140  sm 

1141  ad 

1142  [sm 

1143  ad 

1144  tm 

1145  ad 

1146  ts 

1147  ad 

1150  tm 

1151  ad 

1152  ts 

1153  ad 

1154  tm 

1155  ad 

1156  ts 


35   XJ* 
1 

30   1251 
35   1251 


1127  ad   35   1237 


36  qJ 
1 

35  1241 

25  X  J 

30  2007 

30  2000 

35  2000 

25  2006 

35  2001 

30  1111 

35  1276 

30  2002 

35  1177 

34  1300 

35  2100 
17  2010 
35  1250 

34  1176 

35  2100 
17  2020 
35  1163 

34  1302 

35  2100 
17  2040 


2007 
1235 
2003 

1273 

•  . 

J 

1 

1252 

GJ 
1274 
1274 
2007 
1275 
2006 
1242 
2003 
1277 
1251 
1251 
1177 
3000 
1175 
1250 
2100 
3000 
1247 
1247 
2100 


2003 

1273 

1130 

2007J 

1132 

2000J 

2007 

2000 

2001 

2006 

2002 

2003 

1142 

1251 
1107 
1251 
1164 
1250 
1160 
1250 
1166 
1247 
1160 
1247 
1170 


X^  X  AX^  -  X  3  "/  1  __*.   X J 

i  X  2-24  from  i  X  2~12 

i  X  r12  ■/    i  X  2-24 

£qJ]     in  Mx  and  [*/j  in  M2  of  1130 

[X/  j       in  Mi  and  [dj  in  M2  of  1132 

i\  x  d   x  s"1 


X8 


-1 


3  ft       X  8"*1 


(V'J        in  H3  of  2003 
[cjj'J       in  M3  of  1142 

(i/l)  X  2"12  -*  i  X  2"12 
test    (n^l)  >  i 


1X12" 


12 


i  X  2 


,-12 


print  j  X  2"24  if  ts  2010  up 
(i/l)  X  2"24~j  X  2-24 
test  4  >  j 
0  ->  j  X  2"24 
accumulate    X    if  2020  up 
k-/  1  -+  k 
test  K  >k 
0  -*k 

print  X±  octal  if  2040  up 
53 


1157 

tm 

34 

3000 

2100 

1050 

U.  T,  to  new  print  routine 

1160 

U.  T.  to  1 

L  SR 

1161 

1162 

1163 

f 

f 

0001 

1  X  2"36 

1164 

pr 

21 

1250 

2100 

0001 

print  j  X  ; 

r24 

1165 

tm 

34 

3000 

2100 

1147 

U.  T. 

1166 

ad 

35 

1301 

1304 

1304 

X0/  (m/Df— -  V/  nf  =  t 

1167 

tm 

34 

3000 

2100 

1153 

U.  T. 

1170 

ts 

17 

2020 

3000 

1173 

1171 

(> 

21 

1305 

2100 

»1 

print  Xx,  X2;  . ...  ;Xn 

1172 

tm 

34 

3000 

2100 

1160 

1173 

[pr 

21 

1304 

2100 

n/l] 

print  XQ  - 

t,  X-j^,  Xg,  . . .  «xn 

1174 

tm 

34 

3000 

2100 

1160 

1175 

f 

0001 

f 

1  X  2-24 

1176 

f 

0004 

f 

4  X  2"24 

1177 

0001 

f 

f 

1  X  2'12 

1200 

ad 

35 

1230 

1300 

2000 

f*.]  ** 

Mi  of  2000 

12C1 

ad 

35 

1231 

2000 

1232 

IKl   * 

Hi  of  nod.  dummy  1 

1202 

ad 

35 

1300 

2000 

2002 

[<U   ta 

Mj  of  2002 

1203 

ad 

35 

1233 

2002 

1234 

ft  1    in 

IL  of  mod,  dummy  2 

1204 

sm 

30 

2000 

1235 

2003 

[*0    ** 

LI2  of  2003 

1205 

ad 

35 

2002 

2003 

2004 

[%]  in 

Mi  and  [X  J  in  Mg 

1206 

ad 

35 

1236 

2004 

1237 

[%1  * 

M-,  and  [**]in  1L, 

1207 

ad 

35 

1240 

2000 

1241 

U]  ta 

M]_  of  nod.  dummy  4 

1210 

sm 

30 

2002 

1242 

2002 

M  in 

M,  of  2002 

2004 


54 


1211 

sin 

30 

1300 

1242 

2000 

n  ■/■  lX2"ib  in  2000 

1212 

ad 

35 

1243 

2000 

1173 

sots  up  print  of  XQ#  X-^,  .... 

1213 

su 

36 

2000 

1163 

2001 

n  V(  2"36  in  2001 

1214 

ad 

35 

1244 

2001 

1171 

sets  up  print  of  Xn,Xg  ,  .... 

1215 

ad 

35 

2100 

2100 

1270 

0  -»  i  x  2"36 

1216 

ad 

35 

2002 

1270 

2005 

[^J     in  II3  of  2005 

1217 

ad 

35 

1245 

2005 

1220 

[c\,~]      in  M3  of  1220 

1220 

[ad 

35 

2100 

2100 

*] 

o-qi 

1221 

ad 

35 

1163 

1270 

1270 

(i/l)  X  2-36  _^i  x  2"36 

1222 

"bin 

34 

1270 

2000 

1224 

test  i  >  n  /l 

1223 

tm 

34 

3000 

2100 

1216 

1224 

ad 

35 

2100 

1163 

1303 

signal  for  first  time  through 

1225 

tm 

34 

3000 

2100 

1104 

U.  T. 

1226 

pr 

21 

2006 

1227 

0001 

print  r^ 

1227 

ht 

22 

f 

f 

f. 

halt 

1230 

1304 

f 

f 

W  X  2-12 

1231 

su 

36 

f 

2001 

2004 

dummy  1 

1232 

[su 

36 

*o 

2001 

2004~j 

modified  dummy  1  for  1112 

1233 

sia 

30 

f 

1253 

2001 

dummy  2 

1234 

[sm 

30 

% 

1253 

2001^] 

modified  dummy  2  for  1111 

1235 

02 

f 

f 

0014 

1236 

su 

36 

f 

£ 

2007 

dummy  3 

1237 

f~su 

36 

q0 

*o 

2007  J 

modified  dummy  3  for  1130 

1240 

mr 

25 

f 

1263 

2000 

dummy  4 

1241 
1242 

[mr 

25 

02 

• 
f 

1263 
f 

2000J 
0030 

modified  dummy  4 

\ 
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1243  pr 

1244  pr 

1245  ad 

1246  tm 
1247 
1250 
1251 
1252 
1253 
1254 
1255 
1256 
1257 
1260 
1261 
1262 
1263 
1264 
1265 
1266 

1267  mr 
1270 
1271 

1272  ad 
1273 
1274 
1275 


21  1304 

21  1305 

35  2100 

34  3000 


£ 

f 
f 
f 
4000 
5520 
02   5520 
6525 
4000 
5520 
02   5520 
4000 
25   2004 


2100  f 

2100  f 

2100  f 

2100  f 


f 
f 
f 
f 

f 


0001 

f 
f 

0001 

f 


2357  4151 

2357  4151 

2525  2523 

f  f 

2357  4151 

2357  4151 

f  f 

1257  2005 


35   1304   2007   1304 


02     f 
3000 


f    0003 
f     f 


dummy  print  for  1173 

dummy  print  for  1171 

dummy  for  1220 

dummy  exit  command 

k  X  2"36  (print  tally) 

j  X  2"24   (iteration  tally) 

i  X  2"12  (var.  no.  tally) 

i  X  2"12  -f    j  X  2"24 

B° 

B1 

B2 

B3 

A0  =  1/2 

A1*  /V? 

A2  =  xto 

A3  =  5/6 
C°  =  1/2 

c1  *  m 

c2=  -%/T 
c3  =  1/2 

dummy  for  1115 

i  X  2'36 

i  X  2"12  ■/  i  X  r36 

dummy  for  1124 

i  X  2"12  /  i  X  2"24 

3  X  8""1 
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1276  sm   30   2002   1277     f 

1277  f     f    0003 

(Cells  1300  through  1332  were  filled  by  tape*  These  cells  contain  the 
values  of  t,  z^,  w^#  etc  These  values  arc  not  inoluded  here  because  they 
depend  upon  the  values  of  computed  parameters— see  Chapter  V.) 
Some  notes  on  the  program: 

All  addresses  appearing  in  brackets  are  not  put  in  by  punch  cards  but 
rather  by  other  parts  of  the  program.  They  are  included  for  purposes  of 
clarity  in  deciphering  the  program* 
The  following  notes  apply  to  the  integration  sub-routine: 

ts  2010  —  up  for  iteration  number  tally 

ts  2020  —  up  for  generation  of  xQ 

ts  2040  —  up  for  print  of  x^  octal 

n  ■  no,  of  dependent  variables  (equations) 
X  s  incremental  change  of  independent  variable 

K  =  readout  at  every  K  point 

x  ■  independent  variable 

Hew  print  or  readout  routine  to  start  at  1050. 

Automatic  print  of  cell  2006  indicates  interval  of  integration  too 
large. 

Error  of  integration  is  proportional  to  the  5th  power  of  X» 
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CIIAPTKR  V 

CQNFLUDBIG  RELIARKS  AND  RECOMMENDATIONS 
FOR  FURTHER  WORK 

As  noted  previously,  it  would  be  fatuous  at  this  point  to  draw  any 

definite  oonolusions  anent  the  solution  to  our  problem  proposed  in  the 

first  chapter.  On  the  basis  of  the  computer  solutions  we  cannot  say  with 

absolute  certainty  that  the  results  are  compatible  with  experimental 

observations,  YJe  can  say,  however,  on  the  basis  of  our  sketchy  results 

that  the  theory  and  method  are  promising.  Without  equivocation  we  can 

say  that  the  computer  program  is  eminently  successful— it  does  work  and  it 

docs  solve  the  equations  proposed  as  a  description  of  the  traveling 

density  waves.  A  great  deal  of  work  is  forecast  in  order  to  establish 

the  identity  of  these  equations  with  the  phenomenon  in  the  positive  column. 

Great  numbers  of  initial  conditions  should  be  inserted  in  the  program  and 

the  coefficients  should  be  varied  in  order  to  eventually  establish  the 

validity  of  the  equations.   It  is  believed  that  if  this  is  done#  the 

resulting  solutions  by  the  oomputer  will  correspond  to  the  observations 

made  under  actual  conditions  in  the  laboratory.  However,  no  precise 

conclusions  can  be  drawn  from  a  very  few  computations.  These  computations 

cannot  be  made  in  a  few  hours.   It  is  predicted  that  the  collection  of 

enough  data  will  require  on  the  order  of  several  hundred  hours  of  computer 

time.  The  important  point  to  be  remembered  is  that  a  workable  program  is 

available  and  patience  and  perseverance  will  reward  the  eventual  researcher 

with  the  answer.   It  appears  from  initial  computer  runs  that  the  term 

ztiq   is  dominated  in  the  equations  (4)  by  other  terms  when  X>  s  nQ. 

The  initial  oomputer  run's  suggest  that  "floating  point  arithmetic"  might 
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be  necessary  throughout  the  program  due  to  the  drastic  changes  noted  in 
derivatives   and  the  time  required  for  computation  using  the  present 
program.     Ho  other  trends  wore  noted  from  the   initial  o computer  runs. 
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